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1 Introduction 


The anharmonic oscillator H = uja)a + 7 W 4 is a basic quantum mechanical prob¬ 
lem with one particularly interesting feature: its perturbation series diverges, but 
nevertheless there exist eigenstates and energies which are smooth as the (positive) 
coupling constant 7 goes to zero [[id], [ 12 ], [H|. A similar phenomenon is expected to 
occur in many interacting quantum field theories. The anharmonic oscillator can 
in fact be considered as a (0 + 1 )-dimensional <p A ” field” theory with one degree of 
freedom. 

I 11 this paper, we study the analog of this model in the framework of g-deformed 
quantum mechanics, based on the g-deformed Heisenberg algebra introduced in [[|. 
In particular, one would like to know how the perturbation theory of the g-deformed 
anharmonic oscillator behaves compared to the undeformed case. This is of interest in 
view of a possible g-deformation of field theory, which is expected to be less singular 
than field theory based on ordinary manifolds, since g-deformation generically puts 
physics on a g-lattice !,!• With this motivation, we study the perturbation theory of 
the anharmonic oscillator in terms of the g-deformed harmonic oscillator, which was 
introduced in [ 0 , |[] and realized in the framework of g-deformed quantum mechanics 
in [ 0 . 

There is considerable freedom in defining a g-deformed anharmonic oscillator for 
g ^ 1. Taking advantage of this freedom, we show that for a suitable definition 
of the anharmonic oscillator, the perturbation series converges to exact eigenvalues 
and eigenstates for 1 < q < 1.06 with a certain radius of convergence in 7 . In the 
limit q —» 1 , the model reduces to the usual anharmonic oscillator, and the radius of 
convergence goes to zero. The upper limit on q is not significant. 

This paper is organized as follows: In section 2 we review the g-deformed harmonic 
oscillator and its spectrum, and calculate the relevant matrix elements. In section 3, 
the perturbation series for eigenvalues and eigenstates is discussed. Some estimates 
for the matrix elements are given in the Appendix. 
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2 The q-deformed harmonic oscillator 


In this section, we give a brief review of the g-deformed harmonic oscillator, and 
its realization in terms of a g-deformed Heisenberg algebra. For a more detailed 
discussion, see |T| and ||. 

The q-deformed Heisenberg algebra is the star-algebra generated by X, P , U with 
the relations || 

(pXP - q-^PX = iU (1) 

UX = q~ 1 XU, UP = qPU. 

We assume q > 1 to be real. The star structure is such that X and P are hermitian, 
and U is unitary: 

X = X\ P = P\ = U~\ (2) 

This algebra has the following (momentum-space) representation j5|: 

P\n, cr) = aq n \n , a) 

U\n, a) = \n — 1, a) 

f/ -1 |n, a) = |n + 1, a) 

q~ n i _i 

X\n, cr) = ia -— (g 2 \n — 1, a) — q 2 n + 1, a)) 

q-q 

CT ) (3) 

with n,m G IN and cr, cr' = ±1. The completion of these states defines a Hilbert space 

n. 

The two values of a describe positive respectively negative momenta. (j3|) is a 
star-representation, i.e. the star is implemented as the adjoint of an operator, and 
both A" and P have selfadjoint extensions. That is the reason for introducing cr, see 


This is a starting point for studying g-deformed quantum mechanics Bill 
In particular, one can define g-deformed creation and anihilation operators as follows: 


a = aU~ 2M + /3U~ m P (4) 

a t = aU 2M + fiPU M 
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with M G IV, and a, (3 G C. They satisfy the Biedenharn-Macfarlane algebra [|2j. |3|: 

aa t — q~ 2M a l a = (1 — q~ 2M )aa = 1 (5) 

where we fix a = —--. The occupation number operator is dehned as 

yj\- q ~ 2M * 

h = a f a = aa + (3fiP 2 + a(i{U M + q M U~ M )P. (6) 

Now one can write down the following Hamiltonian, which constitutes the g-deformed 
harmonic oscillator: 

H 0 = ua^a (7) 

The spectrum of H 0 acting on H consists of a bounded spectrum with eigenvalues 
E= uj [n] m = oj 1 1 ~^- 2 M which is 2 M-fold degenerate, and an unbounded spec¬ 
trum with eigenvalues uj(q 2mM + \Z q q - 2 m )• The 2 M ground states of the bounded 
spectrum are 

Q° / \ 71 

|0)‘">= E -4 9+ 

n=—oo V P J 

0 < n < M. (8) 


The existence of an unbounded spectrum beyond is clear in view of (|^), 

since P is an unbounded operator on Pi. For simplicity, we will only consider M = 1 
from now on, and omit the labels /i and M. 

So far, (3 was arbitrary. Requiring that the a, a ^ are smooth for q —> 1 and become 
the usual (undeformed) creation and anihilation operators in the limit, one finds |lj] 


that 


\/l — q 2 ’ ^ y/2muj 
where m is the mass. For this choice, H 0 can be interpreted as a g-deformation 
of the usual harmonic oscillator, and this will be understood in the following. The 
normalized states of the bounded spectrum are 


n) a = 



( 10 ) 


where [n] = ■ We define H b,± C 7i to be the closure of the space spanned by 

the |n)±i. As q —> 1, 7i b,+ becomes the Hilbert space of the usual harmonic oscillator, 
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while the unbounded spectrum disappears at infinity, and the support of the states 
with a — —l goes to —oo in the momentum representation. We will thus concentrate 
on H b,+ . 


The eigenstates of Hq can also be written in terms of the g-deformed Hcrmite 


polynomials, which satisfy (see |10 


fit?® = + 2 9 - 2 [n]/fS(f)) 


( 11 ) 


Defining £ = yjmuX , one has 

\n)c = -4— 

V 2 [nj! 

Using these Hcrmite polynomials, it is straightforward to calculate the action of 
X on an eigenstate \n) a , and it follows in particular that X ■ 7i b,+ C 7i b,+ . This will 
be important for the perturbation theory below. 

Now we turn to the anharmonic oscillator. The undeformed anharmonic oscillator 
is defined by H = cucda + 7 A " 4 for 7 > 0, thus one might naively take the same 
expression for q > 1, and study its perturbation theory. The relevant matrix elements 
can be calculated e.g. using ([[]]), and we find the following results 0 : 

(n|A 4 |n) = ( ( ? 8 ” +6 (J n + 1 K[ n + 2] + g” 4 [n + 1] + g~ 8 [n]) 

+ q~ 8 [ n \([ n + 1 ] + - 1 ])) 

(n + 4|A" 4 |?i) = ^ (f n+1A \] [n + 1] [n + 2] [n + 3] [n + 4] 

<n + 2|A 4 | n) = q 8n+12 \/[n + l}[n + 2] 

([n + 3] + q~ A [n + 2] + q~ 8 [n + 1] + q~ 12 [n]) (12) 


They are independent of a which is suppressed. All other nonvanishing matrix ele¬ 
ments can be obtained from those by hermiticity. 

Looking at the powers of q in the matrix elements, one quickly finds that the 
perturbation series diverges even faster than in the undeformed case. 

However, it is important to realize that there is no reason for considering the same 
expression for H as in the undeformed case; the only requirement one has to impose 
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is that H should reduce to the usual anharmonic oscillator as q —» 1. Therefore we 
might just as well consider the Hamiltonian 


H = H 0 + 7 H' 


(13) 


with 

H'= ^(X 4 Q 5 + Q 5 X 4 ), where 
Q = (1 - a f a(l - q~ 2 )). (14) 

Q satisfies 

Q\n) = q~ 2n \n). (15) 


The matrix elements (n\H'\m) can be easily obtained from (p~2D , see Figure 2. As is 
shown in the Appendix, they have the following upper bound: 


(n\H'\m) < C(q) 


[3] 4 [2] 8 ,- 2 ”“ +1 > ma J 2 


V°[3] 4 [2] s 
81(1 — <?- 2 ) 2 


(16) 


for 1 < q < 1.06, where n max = In view of the results of the next section, we 

define (|T3|) to be the g-deformed anharmonic oscillator. 


3 Perturbation Expansion 

We will use the standard Rayleigh-Schrddinger perturbation formulas for the eigen¬ 
states and eigenvalues of 


H = H 0 + H x = H 0 + 7 H’ 
in terms of the unperturbed ones, H 0 \n) = E^\n): 


(17) 


OO OO / 1 

A E n = Y, E { n k \AE n , 7) := ( Z7nT ~ Qn (^i - A E n ) ) | n) 

k=0 k=0 


Pl(0) TT 

. fin — n 0 


(18) 
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where Q n = (1 — \n)(n\), and 


I E n ) = 


In 


oo K> n - AEn\nj) (n k \Hi\n) 

> + E £ —^; ' 


k=l 

nr^n 


n (^ 0) - Enj 

3 = 1 


( 0 )'. 

i J 


(19) 


Strictly speaking, we are of course dealing with a degenerate problem (since a = 
± 1 ); however as already explained, A" and Q leave 7i b,+ invariant, thus the two values 
of a do not interfere, and we can restrict ourselves to the a = +1 sector. This will be 
understood in the following. We will show that these series in fact converge to exact 
eigenvalues and eigenstates of the q-deformed anharmonic oscillator, for a certain 
range of 7 which depends on q. 


3.1 Energy Levels 


If 7 and A E n are not real, then Hi is understood to act on the right in the above 
formulas, so that the matrix elements can be continued analytically in 7 and A E n . 
We show first that the sum in fll 8 |) is absolutely convergent for \AE n \ < u;/5 and 
| 7 1 < 7 (q), where 7 (q) > 0 provided q > 1, see (j 22 |). Thus the rhs of ( |T 8 | ) is an 
analytic function of A E n and 7 in that domain, which can be solved for A E n by the 
implicit function theorem, defining an analytic function A E n {p/). 

To see that the sum in (|1^) is (absolutely) convergent for a certain range of A E n 
and 7, we hrst write E^ more explicitely: 



bW(AB„, 7 ) 


(n\Hi\n) 

^2 (n\ H i\ni){ni\ H i\n) 


"1 (Ei 


( 0 ) 


ndo) i 

■Eni . 


<k -1 


(n\Hi\ni) n {n j - 1 \H 1 - AE n \nj) 

E ; 


^—' / rdO) ztHuj \/z7H u j 7 tA u J \ ( nA u ) cH u A 

71,»2,...,n fe _ 1 ( h, n — — lhn k _ 2 ) . . . \Ihn ~ Hini ) 

n r ^n 

for k > 3 


do) 


do) 


do) 


do) 


dob 


( 20 ) 
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As is shown in Appendix A, the following estimate is valid for q e]l; 1.06[: 

I^v)i < == 

for k > 2 , 


( 21 ) 


The factor 5 comes from the fact that for any given rij , there are only 5 possible 
rij + 1 such that the matrix elements in the perturbation expansion do not vanish (see 

my 


The series fllq) is absolutely convergent if the following condition holds: 

Ei k+ V 


n(k) 


< 6 for some 6 < 1 


Now 


J 7 |C(?) + |AB„| ^ 

= 5-—--- < 0, 


i(fc) 


[2] eng 


—2 n 


and we find that the condition holds e.g. for \AE n \ < tn/5 and 

tn([ 2 ]q~ 2n - 1 ) 


M < 7(g) : = 


5 C(q) 


( 22 ) 


Therefore we have shown that in this domain, the rhs of (p~8|) dehnes an analytic 
function in A E n and 7 . Notice that 7 (g) —> 0 as g —> 1. 

Now consider the equation 

OO 

G(AE n , 7) := £ E£\AE n , 7) - A E n = 0. 

fc =0 

In the above domain, this is a uniformly convergent series of analytic functions (for 
fixed q in the interval ] 1,1.06[, say). But then using (PU|), one sees that 


dAE, 


E£o EP(AE n , 7 ) 4 _. o = 0 , i.e. 


7=0 


d 


dAE, 


-G(AE n , 7) 7 ^ 0 


(23) 


for 7 and A E n in a neighborhood of 0, by analyticity. Now the implicit function theo¬ 
rem states that there is a function AE n {p/) which solves (|25|) and satisfies AE n {{)) = 0. 
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Moreover, AE n {^) is analytic in a neighborhood of 0, and \AE n \ < u;/5 holds auto¬ 
matically if 7 is small enough. 

The domain of convergence 7 (g) is shown in Figure 1 for q e]l; 1.06 [ and uj — 1. In 
particular, 7 (g) goes to zero for q —> 1 , in accordance with the well-known fact that 
the perturbation series for the undeformed anharmonic oscillator is divergent |T 2 . 


3.2 Eigenstates 

In this section, we show that (|19|) converges in H b,+ C Tt for I 7 I < 7 (g) and 1 < 
g < 1.06, where A E n = AE n {^/) is now the perturbed energy found in the previous 
section. To do this, we have to show that 

OO 

E IHK)| 2 <oo, (24) 

m =0 

or more explicitely 


(E n \E n ) = £ \{m\E n )f 


m =0 


= E 


m =0 


<5 m ,n + E E 


5m,ni In (n j -i\H 1 - AE n \nj) \ (n fe |i7i|n) , 2 


k=l n i’- - n k 

rif^n 


n {e£ ] - e^) 

j =1 


From the form of the matrix elements (0), we see that the second term is nonzero 
only for k > therefore 


(E n \E n ) < 1 + E E 

m=0 1.^ \m- 


k>l . 

— 4 


<i + E h 


m =0 


^\C(q) + \AE n \) k -^ \C(q)5 
([ 2 ]g- 2 "w) fc 

\m—n\ 

~ 'E(s 

vk =0 


rfc\2 


h\c(g) + \AE t 

[2}q- 2n co 


| 7 | C(q) + |Aff„ 

[ 2 ]g _ 2 n u; 



for 1 < g < 1.06. Clearly this converges for 7 in the analyticity domain defined 
above (such that \AE n \ < u/h as before), therefore the series (pT|) converges in 7 i b,+ . 
Finally, both El 0 and H' leave 7i b,+ C Et invariant and are bounded operators on 
H b,+ (H' is bounded because of (|27|) and the fact that H' acting on | n) has no more 
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that 5 nonvanishing components in terms of that basis). Now it follows that \E n ) and 
E® + A E n are indeed eigenstates and eigenvalues of the full anharmonic oscillator. 

As already mentioned, it is known [|T^| that the undeformed anharmonic oscillator 
does have nonperturbative eigenstates and energies for 7 > 0 , which are nevertheless 
smooth as 7 goes to zero from above. Now the formulas (|T§| ) ff. can be analytically 
continued in q as well, and one would expect that the above domain of analyticity for 
A E n and 7 can be extended to include q = 1 and positive real axis of 7 . However, at 
present we are not able to show this. 

Acknowledgements A. P.N. acknowledges with thanks the support from MPI. 


Appendix: Matrix elements 

Because [n] is an increasing function in n, we have the following estimates: 


k(l + (f 40 )9 14_2 >] 2 < (n + 4| H'\n) < b(l + g - 40 ), 14 - 2 > + 4 ] 2 
b(l + ?- 20 ) 9 12 - 2 ”[ 4 ] 4 [r ,] 2 < (n + 2| H'\n) < ]p(l + ?' 2 < >)<( 12 - 2 ”[4] 1 [n + 3 ] 2 


7 V 2 “ +e [ 3 ] 4 [ 2 ] s [n ] 2 < (n| H'\n) < 7 9 - 2 ” + «[ 3 ] 4 [ 2 ] 8 [n + 2 ] 2 

(25) 

with 

1 -q~ ni 1 -q~ 2n 


See Figure 2 for a plot of (n\H'\n). 

To simplify this, consider the function q~ 2n [n] 2 for n E JR, which takes its maxi¬ 
mum value 81 ( 1 _ 4 g -2 )2 at ra = n max , 


Tlmax 


In 3 
2 In q 


(26) 
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Figure 2: The matrix elements (n\H'\n) for q E [1.001,1.002] depending on n 

The matrix elements have a maximum for n close to n max . More precisely, we can 
show the following estimate: 

I (n + i\H’\n)\ < C(q) := ,-2«~+ 1 »[ 3 ] 4 [2] 8 [n mM! ] 2 = g^^ 13 ^ 12 ]^ (27) 

for all n, m E IN. Indeed, 

<n + 4|if» < ^(1 + q~ A0 )q 14 ~ 2n [n + 4] 2 

= i(l + g- 40 )g 2 V 2(n+4) [n + 4] 2 

< ^(1 + q~ 40 )q 22 q- 2nmax [rimax] 2 , 

furthermore 

(n + 2\H'\n) < i(l + q- 20 )q 18 [A} 4 q- 2 ^[n + 3] 2 

< + q- 20 )q 18 [4] 4 q- 2n — [n max ] 2 , (28) 
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and 


(n\H'\n) < 9 W (3] 4 [2] s ? 2n ~[n m „] 2 = C(q ) 


(29) 


Now for 1 < q < 1.06, one has 


2<T 16 [3] 4 [2] 8 
1 + g- 40 


(30) 


(for i — 4) and 


2 [3] 4 [2] 8 

? 12 [4] 4 (l + g- 20 ) 


(31) 


(for i — 1). Combining these estimates, we obtain (|27|) . Furthermore |A!® — e !£±J > 

[i]q~ 2n co, therefore |— E^\ > [2 ]q~ ~ 2n ui in the denominators of the perturbation 
expansion, since i > 2. Now m follows, because for any given rij in the perturbation 
series, there are at most 5 possible n. J+ \ such that (rij\H l \rij + i ) is nonzero; this means 
that the number of terms at order k is at most 5 fc ~h 
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